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/A 


Direct conversion of thermal energy into electrical energy using a photovoltaic cell is called 
thermophotovoltaic energy conversion. One way to make this an efficient process is to have the 
thermal energy source be an efficient selective emitter of radiation. The emission must be near 
the band-gap energy of the photovoltaic cell. One possible method to achieve an efficient 
selective emitter is the use of a thin film of rare-earth oxides. The determination of the efficiency 
of such an emitter requires analysis of the spectral emittance of the thin film including scattering 
and reflectance at the vacuum-film and film-substrate interfaces. Emitter efficiencies (power 
emitted in emission band/total emitted power) in the range 0.35-0.7 are predicted. There is an 
optimum optical depth to obtain maximum efficiency. High emitter efficiencies are attained only 
for low (<0.05) substrate emittance values, both with and without scattering. The low substrate 
emittance required for high efficiency limits the choice of substrate materials to highly reflective 
metals or high-transmission materials such as sapphire. 


I. INTRODUCTION 

Thermophotovoltaic (TPV) energy conversion re- 
search dates back to the 1960s. 1,2 Most of the early 
studies 3-7 of TPV energy conversion used the solar flux as 
the thermal energy input. In those systems the solar energy 
was absorbed and reemitted, the idea being to improve the 
efficiency of the photovoltaic (PV) conversion by match- 
ing the spectrum of the radiation to the band gap of the PV 
cell. There is now renewed interest in TPV energy conver- 
sion for space power systems using nuclear energy 
sources 8,9 and also commercial combustion driven 
systems. 10 

The preferred emitter for the early solar driven 
systems 3-6 was a blackbody coupled with a bandpass filter. 
Ideally, only photons with energy greater than the PV cell 
band-gap energy are transmitted and all others are re- 
flected back to the emitter to be reabsorbed. Thus, the 
emitter-filter combination performs as a selective emitter. 
For such systems to be efficient the filter absorptance must 
be low. 4,6 A simpler, efficient TPV system will be possible 
if an efficient selective emitter can be produced. In Ref. 7 a 
selective emitter is proposed for a solar TPV system. 

The early work of White and Schwartz 11 recognized 
the benefits of selective emitters for efficient thermophoto- 
voltaic energy conversion; however, finding an efficient se- 
lective emitter has been a difficult task. The most promis- 
ing solid selective emitters have been the rare-earth 
elements. 12 For doubly and triply charged ions of these 
elements in crystals the orbits of the valence 4 / electrons, 
which account for emission and absorption, lie inside the 5s 
and 5 p electron orbits. As a result, the rare-earth ions in 
the solid state have radiative characteristics much like they 
would have if they were isolated. They emit in narrow 
bands rather than in a continuum as do most solids. The 5s 
and 5 p electrons “shield” the 4/ valence electrons from the 
surrounding ions in the crystal. The spectra of these rare- 


earth ions in crystals have been extensively studied. Most 
of this work is summarized in the text of Dieke. 13 

Early spectral emittance work 12 on rare-earth oxides 
that are suitable for TPV showed strong emission bands. 
However, the emittance for photon energies below the 
band gap for PV materials was also significant. As a result, 
the efficiency of these emitters was low. In the last few 
years, however, Nelson and Parent 14,15 have reported a 
large improvement in rare-earth oxide emitters. Their 
emitters are constructed of fine (5-10 /xm) rare-earth ox- 
ide fibers similar to the construction of the Welsbach man- 
tle used in gas lanterns. The very small characteristic di- 
mension of these emitters results in a low emittance in the 
low-energy part of the spectrum, thus giving them a high 
efficiency. Chubb 16 has calculated the efficiency for the 
rare-earth emitters constructed with a small characteristic 
dimension and found that efficiencies approaching 80% are 
possible. 

As already mentioned, major improvement in the rare- 
earth oxide emitters resulted from the small characteristic 
dimension possible with the mantle-type geometry. An- 
other way to obtain a small characteristic dimension is to 
use a thin film. This geometry enables a smaller ( < 10 /xm) 
characteristic dimension and thus higher efficiency. Also, a 
thin film on a strong substrate is more rugged than the 
mantle-type emitter. 

To evaluate the potential of a thin-film selective emit- 
ter, the analysis described in this article was carried out. 
An expression for the spectral emittance of a thin film was 
derived using a continuum radiative transfer analysis that 
includes scattering. Using the derived spectral emittance, 
the emitter efficiency and power density were calculated. 
To model the photon energy dependence of the rare-earth 
oxide optical properties, the spectrum was split into three 
regions. Within each region the optical properties are as- 
sumed constant. 
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FIG. 1. Schematic of emitting film. 


II. SPECTRAL EMITTANCE OF THIN FILM 


In order to determine the efficiency of a thin-film se- 
lective emitter, the spectral emittance must be known. The 
following deriviation is based on the continuum radiation 
transfer equations. 17 Figure 1 shows the one-dimensional 
model of the film. The radiation intensity i v is assumed to 
vary only in the x direction. This is an excellent approxi- 
mation for thin films where the film dimensions perpendic- 
ular to the x direction are much larger than the film thick- 
ness d. The radiative transfer equations for the spectral 
radiation intensities in the positive x direction i+ and the 
negative x direction i ~ , as functions of the optical depth 
K v and angle 9 (fi = cos 9) are the following: 

C (A' v ,/i)=C (0,/x)expf -- S V (K* ,/i) 

t K v -Kt\ dfCt 

Xexpl— — — — -j (1) 

i~ (K v ,n)= i~ ( K vd ,n ) exp 



-I 


K vd {Kt-Kt\ did 

S v (Kt .M) ex PI 


— ( 2 ) 

And the source function S v for isotropic scattering satisfies 
the following equation: 

s v (* v )=n}(i-a v )/ v6 +y [ J o C (0,/i)exp^-^j^ 

r 1 _ (K vd -K v \ j 

+ J o iy (K vd ,-[i)ex pi — Jdfj. 


+ 


C vd S v (Kt)E ] (\K^-K v \)dKt 
Jo 


(3) 


where is the index of refraction of the film. Also, ap- 
pearing in Eqs. ( 1 )— (3) is the optical depth K v> 


K v (x)= f (a v -\-a v )dx*— f a v dx * (4) 

Jo Jo 

where a v is the spectral absorption coefficient, o v is the 
scattering coefficient, and a v =a v +a v is the extinction co- 
efficient. Also, the scattering albedo ft v is defined as fol- 
lows: 


ft. 


a v +<7 v 


(5) 


and the blackbody intensity i vb and hemispherical emissive 
power e vb are the following: 


e vb 2hv* 1 

* vb ~ it cl zxp(hv/kT) — \ 9 


( 6 ) 


where h is Planck’s constant, k is Boltzmann’s constant, 
and c 0 is the vacuum speed of light. Finally, the general 
exponential integral E n (u) is the following: 



The following approximations were made to simplify the 
solution to Eqs. (l)-(3); 

(1) The temperature T E is uniform; therefore, a v , a v , 
and rif are constant throughout the film and K v =a^x and 
E vd aji. 

(2) The film is a pure dielectric; therefore n f is a real 
number. 

(3) Scattering is isotropic. 

(4) Boundaries at x = 0 and x=d behave diffusively; 
therefore, i+(0), i~{K vd ), and reflectivity at substrate-film 
interface p VJ , reflectivity at film- vacuum interface p v 0 , and 
emittance of substrate e vs are independent of 9. 

(5) Interference effects are neglected although d<X 
for part of spectrum. Assuming diffuse boundaries yields 
the following boundary conditions at x=0 (K v =0) and 
x=d(K v =K vd ), 


q„ (0) = e v! e vb +p v ji v (0), K v = 0, 

(8) 

q;(K vd )= PM q + (K vd ), K v =K vd , 

(9) 

where 


f-rr/2 

q+ (0) =2ir /+(0)cos 8 sin 6 d6 

Jo 


= 2ir f P i+(0)d P =m+{0), 

Jo 

(10a) 

and, similarly, 


q~(0)=m'~(0), 

(10b) 

q+{K vd )=irit(K vd ), 

(10c) 

q v iK vd ) TTi v (K vd ). 

(10d) 


Since i\t (0) and i v ( K vd ) are assumed independent of 
p, the p integration in Eq. ( 3 ) can be performed, 
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S v (^ v )=n}(l-ft v )— +-f 

J 7T L 


E 2 {K V ) 


gy + (0) 

7T 


+E 1 (K vd -K v ) 


V ( K vd ) 
7 T 


+ 


r 


S v (K*)E x (\K*-K v \)dK*. 


(ID 


This must be solved for S V (AT V ) which can then be used in 
Eqs. ( 1 ) and (2). The method outlined in Ref. 8 was used 
to solve Eq. (11). Replacing the exponential integral func- 
tions in Eq. (11) with exponential approximations allows 
the integral equation to be converted to a linear second- 
order differential equation, which can be easily solved. De- 
tails of the solution are given in the Appendix. 

The hemispherical spectral emittance is defined as fol- 
lows: 




Qv(Kvd) 


Cvb 


( 12 ) 


where Q v (K vd ) is the spectral emissive power at K v =K vd . 
To include the refractive properties of the film-vacuum 
interface the maximum angle 0 M that allows radiation to 
escape the film is given by Snell’s law, 


n 2 M =<x>s 2 e M =\- 



(13) 


For 9>0 M radiation is totally reflected. Including angles 
greater than 0 M would allow the possibility of e v >l. 17 
Therefore, using the boundary condition given by Eq. (9), 
the emissive power becomes 


QyiK vd ) 


= 2tt 

Jo 


[C (AT W ,/X) — / v (K vd ,p) ]cos 6 sin 6 dO 


— 27r( 1 — 



j fii + {K vd ,n)dfi. 


(14) 


Now, using Eq. ( 1 ) and the boundary condition given by 
Eq. (8) in Eq. ( 14), 


Q v {K vd ) = (\-p^) 


<lv (**) 2 [ e v5 e v £+ (0) ] 



The Qy ( K vd ) and q v (0) terms are obtained from Eqs. ( 1 ) 
and (2) and the boundary conditions given by Eqs. (8) 
and (9), 

(Kvd) =2tt f fiij (K vdi fi)dfi 
Jo 

— 2 [e v A'6+p v s£ v (0) \K^(K vd ) +0+ (16) 

q~(0)=2rr f \ii~ {Q,fj,)dfj, 

Jo 

= ( K vd ) E, ( K vd ) + d> _ . (17) 

The source term integrals appearing in Eqs. (15) to (17) 
are 


<D + = 2ir \ Kvd Sy{Kt)E 2 (K vd -K*)dK* t (18) 

Jo 

ft* (Kyd \ 

Sy(Kt)E 2 [— Wj, (19) 

Jo \ fl M ) 

<p_ = 27r [ Kvd Sy(K* y )E 2 (Kt)dK*. (20) 

Jo 

Using Eq. (A8) (see the Appendix) for S v in Eqs. (18)- 
(20) and the exponential approximation for E 2 [Eq. ( A2)] 
results in the following after a great deal of algebra: 

< &M=h(f J 'M)Qv (®) + ^(/ i Af)<?v (Kvd) +h(f l M) e vbJ (21) 

< J > + =/iU)9v + (0)+/ 2 (l)^(A:v rf ) +/ 3 OK*, (22) 

O_=/ 2 (l)^(0)4/,(l)^(A: vrf )4/3(l)e v6 , (23) 

where 


l l {v)=2fi 2 Y^j{A A [m 2 (fJ.)A 2 -m l (iJ,)A l ] 

! 2 (n)=2n 2 y^j{A 4 [m l (^)A 2 -m 2 (n)A { ] 

+ m 0 (n)}, 

W (^)[ 1 — 


vd 


a 2 

— Ay4 4 [m { (fi) +m 2 (n)] — 

°2 


m,(/i) = 


m 2 (n) = 


1 +fi( ^/l— z/u) 


1 +fi( ijl —z/u) 


1 —exp I — r— 

V M 


(e-'-e-W), 


m^n) =j— ( e y —e y/>1 ) for 

1 fl 


m 3 (\)=ye y , 

1 


m 0 (p) = 


1+M 


[-expf-^) 


U) , w_ 

A \=~T7 e T e ~ 

«+l u— 1 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 


w_ 

^ 3 = ( ! _, V j (<"+ -w_e~ r ). 


A *=,J 


w 


= 1+ Vl ~Z, 


(33) 

(34) 

(35) 

(36) 
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w_ = 1 — {\ —z, 

g^h 

r_ 2{l-[z/(l-« 2 )]}’ 
r=b\ yjl — zK vd , 
y=b 2 K vd , 
a ift v 

z =—j— > 


(37) 

(38) 

(39) 

(40) 

(41) 



(42) 


After using the boundary conditions [Eqs. (8) and 
(9)] in Eqs. (22) and (23) and substituting the results in 
Eqs. (16) and (17) we are left with two algebraic equa- 
tions for ( K vd ) and q~ ( 0 ) . Solving these yields 


^ +( ^ )= DiN (6vlLl(1) + /3(1) 


X{l+p w [ii(l)-/ I (l)]». (43) 


q~ (0) =5^ {<U L 2 ( l; )/ 2 ( ! ) + PyC L\{ (1 1 ) ] + / 3 ( ! ) 


X[L,(1) + L 2 (1)]}. (44) 

where 

Z, 1 (/i)=/ 1 (/x)+2 M 2 £ 3 ^J, (45) 

L 2 (/x) = l— P vo/ 2 (m), (46) 

DEN=L 2 (1) [1 — p vi / 2 ( 1 ) ] — pvop vs Li( 1 ). (47) 


Using Eqs. (43), (44), and (21) in Eq. (15) the hemi- 
spherical spectral emittance is obtained, 

Qy(K vd ) 

€v e vb 

1 — PvO 

^-^[^s[L l {l)L 2 (n M )-L ] (^ M )L 2 (l)} 


+ hi 1 + Pvs{ [ L 2 (Pm) — pv 0^1 ip\i) ] 


X[I 1 (1)-/ 2 (1)]-IiW})-/ 3 WDEN]. (48) 

In Eq. (48) the contributions to e v are split into two 
groups. The first group, which is proportional to e vs , re- 
sults from the substrate emittance. The second group, 
which is proportional to / 3 , is the contribution of the film 
to the spectral emittance. 

Consider two limiting solutions to Eq. (48). First, for 
K vd ^0 expand Eq. (48) in terms of K vd and retain only 
the linear K vd term. In this case the following is obtained: 




(1 — Pvq)U~Ma/) 

1 — PvOpYS 


e vsi 1 ) + 



X ( 1 —a 2 z){\ +p vs [ 1 1 + Pvo) ]} 



Optical depth, = (ct v + a v )d 

FIG. 2. Spectral hemispherical emittance e v of emitting film for various 
values of scattering albedo ft v , substrate emittance e VJ — 0.1, and no re- 
flection at interfaces (pvo = Pvj— 0). 


-e vj [2-a 2 n v ( 1 -Pa/Pvo) ] 



(49) 


Second, for K vd — > oo Eq. (48) becomes the following: 


1-PvO 

l-Prf^O) 




X t / ? (/ij#)/ 3 " ( 1 ) ■ - /? ( 1 ) If inu) ] }, (50) 

where 


/« 

l 2 


(fi)=2fi 2 



X 


1- 




W + (l—U)[l+p( yjl—z/u)] 

i-a v 


(51) 


2 / i 




2 z 

X 7 . (52) 

w+(\—u)[l+p( 

Equation (49) shows the importance of the substrate emit- 
tance € vs . For small K vd the film emittance is dominated by 
€ vs . In order to have an efficient thin-film selective emitter 
the substrate spectral emittance outside the desired emit- 
tance band must be small. This is discussed further in the 
following section. For very thick films the substrate emit- 
tance makes no contribution to the film emittance, as Eq. 
(50) illustrates. 

Figure 2 shows e v vs K vd [Eq. (48)] at several values of 
the scattering albedo fl v for a substrate emittance e vs =0.1. 
Also, reflectivities at the interfaces are zero (p vs —p v0 = 0) 
so that n 0 = n f=n s = 1 and fi M = 0. Values for the constants 
used in the exponential integral approximations [Eq. ( A2)] 
were taken from Ref. 18 (a { =b { = 2) and Ref. 17 (a 2 — i 
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(53a) 


b 2 —j). Also, shown in Fig. 2 is the linear approximation 
[Eq. (49)] f° r ^e case fl v — 0. Obviously the linear approx- 
imation applies only very close to K vd =0. In the emitter 
efficiency calculations to follow the optical depths of inter- 
est are greater than 0.5; therefore, the linear approxima- 
tions cannot be used. For the case H v =l (a v =0), when 
the film emits no radiation only the substrate contributes to 
e v . 

III. EFFICIENCY OF THIN-FILM SELECTIVE EMITTER 

Knowing e v allows the calculation of the efficiency of a 
thin-film selective emitter. The work of Nelson and 
Parent 14 ’ 15 shows that the emission spectra of the rare- 
earth oxides neodymia (Nd 2 0 3 ), holmia (Ho 2 0 3 ), erbia 
(Er 2 0 3 ), and ytterbia (Yb 2 0 3 ) are dominated by a single, 
narrow strong emission band. Therefore, to simplify the 
analysis, the emittance of the film is split into three regions 
of constant optical properties. For the region of photon 
energy {E— hv) below the emission band the emittance is 
£/, which will be a function of the optical depth K h scat- 
tering albedo fl/, substrate emissivity e sh reflectivity of 
film-vacuum interface p ob and reflectivity of film-substrate 
interface p sh 


€i—€i{K h £ij f € sh pQ h p s i). 

Similarly the emittance in the emission band is e g and the 
emittance for photon energies greater than the emission 
band is e u , 

€g=€g{K gt £lg,€ S g f pog,p S g), (53b) 

€u~€u(K U t^lu^su*POu*Psu) ‘ (53c) 

In Eq. (53) K g is the optical depth for the emission band, 
is the optical depth for photon energies greater than the 
emission band photon energy (E> E g =hv g ). Similarly, il g 
and O u , e sg and e su , p 0g and p 0u , as well as p sg and p su are 
the scattering albedos, substrate emittances, film- vacuum 
reflectivity, and film-substrate reflectivity for the emission- 
band and above-emission-band regions, respectively. 

Now define the emitter efficiency as follows: 

emitted power from emission band 

J] E = — - 7 — . (54) 

total emitted power 

Using the three emittances just discussed, the efficiency 
becomes 


J 


gJ E a 


*S+ A V 2 




^Eg/i e Eb 


dE 


e /f 0 8 * Eg/2e E b ^ tfy2 e b:b dE+*uf E g +bEJ2 e Eb ^E 



(55a) 


(55b) 


where e Eb (E>T E ) is the blackbody emissive power [Eq. 
(6)] in terms of photon energy E—hv and emitter temper- 
ature T E . The photon energy at the center of the emission 
band is E g and AE g is the width of the emission band. Also, 
using Eq. (6) for e Eb in Eq. (55) yields 


G(s,t)-- 


s(l+//2) 




{U. 
(/. 

hi; 

u 


s(\-t/l) 


s( 1 —t/2) 


+ t/ 2) 


**-1 


e x — 1 


dx 


dx 


e*-l 




dx 


s(\ + t/2) X 
s(i-t/2) e*— 1 


dx 


(56) 


(57) 


To simplify the presentation of results, define the op- 
tical properties below and above the emission band in 
terms of the emission-band properties. Therefore, for the 


extinction coefficients (a = a-b<r), 

ai=ai+(Ji=fi(ag+ag) or Ki=ffK g , (58a) 

a u —o u +a u =f u (a g +a g ) or K u —fji g , (58b) 

and for the scattering albedo [Cl — a/(a~\-a)], 

tti=gfi g > (59a) 

(59b) 

and for the substrate emittance (e 5 ), 

Csi= hi^sg » ( 60a ) 

€ su h u € sg . ( 60b ) 


As Eq. (55) indicates rf E is a function of Eg/kT E , AEg/E g , 
and the emittance ratios e/€ g and e u /e g . And as Eqs. 
(53a)- (53c) show, the emittance ratios are functions of 
the film and substrate optical properties. For given values 
of E g /kT Eg and A Eg/E g there will be an optimum film 
thickness that will maximize rf E . 


Now the emittance ratios €/€ g and e u /e g that appear in 7] f 
are determined by the emission-band optical properties, the 
interface reflectivities, and the /, g, and h parameters. To 
obtain the maximum efficiency t j e for given values of 

E g /kT E and A Eg/E g and the /, g, and h parameters, the 
optical depth K g must be varied until Eq. (55) is a maxi- 
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mum. This optimization process was carried out on a per- 
sonal computer using mathematica 19 software. 

Consider the case where the substrate, film, and vac- 
uum surroundings all have the same index of refraction 
(n s =rif=n 0 = 1). In this case and p v o=Pvs = ®' The 

optimum efficiency 77 E as a function of Eg/kT E for sev- 
eral values of the substrate emittance e sg are shown for this 
case in Fig. 3. Figure 3(a) is for no scattering (0=0) and 
Fig. 3(b) is for equal scattering and an absorption coeffi- 
cient (fl of 0,5). It is assumed that LEg/E g —0A y which is 
typical for the rare-earth oxides considered by Nelson and 
Parent 14 ’ 15 , gi=g u =ht=h u = 1 and that fi=f u — 0.01, 
since /, g, and h are the same for both above- and below- 
emission-band regions, € u = € h and there is a single opti- 
mum optical depth independent of Eg/kT Ey Kg opt » for max- 
imum 7 ] e , If fjgy and h are different above and below the 
emission band then there will be a different K gopt that yields 
j] E for each value of Eg/kT E . The optimum optical depth 
for each value of e sg is listed in Fig. 3. For Fig. 4 the same 
conditions were used as in Fig. 3 except that /„=//= 0.1. 
In other words, the extinction coefficient below and above 
the emission band is a factor of 10 higher in Fig. 4. 

From the work of Nelson and Parent 14,15 with the rare- 
earth oxides already mentioned, it is expected that f u and 
// should be in the range of 0.01-0.1. Thus, Figs. 3 and 4 
should be representative of the range of emitter efficiencies 
that can be expected from a thin-film rare-earth oxide se- 
lective emitter. For /„= f 0.01 maximum efficiencies up 
to 0.7 are expected and for f u =f e = 0.1 maximum efficien- 
cies of 0.35 are expected. There are several important con- 
clusions to be made based on the results shown in Figs. 3 
and 4. First of all, by comparing Fig. 3(a) with Fig. 3(b) 
and Fig. 4(a) with Fig. 4(b) it can be seen that, for 
0.5, scattering has a negligible effect on the 
magnitude of the optimum efficiency; however, to obtain 
t] e with scattering requires larger optical depths. With 
scattering the same at all photon energies, we expect scat- 
tering to have negligible effect on t} e since it is the emis- 
sive properties (absorption coefficient) that should deter- 
mine t) e . Scattering will reduce the emission (Fig. 2) but it 
will not greatly alter the emittance ratios e u /e g and e/e g 
that determine rj E [Eq. (55)]. Second, the importance of a 
low-emittance substrate in order to have high efficiency is 
evident from Figs. 3 and 4. The range of Eg/kT E for at- 
taining large r] E also decreases for smaller € sg . Also, the 
optimum optical depth increases as e sg increases. As f u and 
fi increase, the effect of the substrate emittance is reduced. 
Finally, the maximum efficiency in all cases occurs at 
Eg/kT E ~\. This result was obtained by Chubb 16 using Eq. 
(55) for a selective emitter with assumed values for € u /e g 
and 6//e g in the range of 0.01-0.2. If f t > f u the optimum 
value for Eg/kT E is reduced (higher T E required for given 
E g ). Table I shows the optimum temperature 
(Eg/kT E — 4) for four of the rare-earth oxides when 
f u =fj. For Yb 2 0 3 and Er 2 0 3 , reaching the optimum tem- 
perature ( T e > 3000 K) means that Yb 2 0 3 and Er 2 0 3 
will melt; however, for Ho 2 0 3 and Nd 2 0 3 the optimum 
temperature should be attainable. 


€ sg* Kgopt 



Eg/kT E 


e sg» Kgopt 



Eg/kTe 


FIG. 3. Dependence of optimum emitter efficiency VE opX > on dimension- 
less emission-band energy, using the optimum optical depth Kg ^ , for 
both below-emission-band and above-emission-band optical depths equal 
to 0.01 (emission-band optical depth) (//=/ u =0.01). Also, substrate 
emittance is constant (h u =h { = 1), no reflectance at interfaces 
{p^ = p vs = 0), and A£y/£y=0.1. (a) No scattering (n / =fl u = n^=0); 
(b) equal scattering and absorption coefficients (n/=ft u =n ? =0.5). 

Figure 5 shows the dependence of efficiency on optical 
depth for no scattering (n/=fl u — fl g =0). In Fig. 5(a) the 
efficiency is shown as a function of K g for the same condi- 
tions as in Fig. 3(a) but with Eg/kT E fixed at the value for 
optimum efficiency (Eg/kT E =4). Figure 5(b) shows t} e 
vs K g for the same conditions as Fig. 4(a), again with 
Eg/kT e =A. When the emission-band extinction coefficient 
is much greater than the below-band and above-band ex- 
tinction coefficients (f u =f t = 0.01) then near-optimum ef- 
ficiency is attained as long as K g > 1. However, for 
f u = f /=0. 1 optimum efficiency is achieved only when 
K g ^\. The optimum optical depths for / w =//— 0.1 [Fig. 
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€ sq> Kgopt 



Eg/kT E 


e sg> Kgopt 



Eg/kT E 


FIG. 4. Dependence of optimum emitter efficiency 7 j E , on dimension- 
less emission -band energy E/kT E , using the optimum optical depth 
K g ( for both below-emission-band and above-emission-band optical 
depths equal to 0.1 (emission-band optical depth) (//=/ M =0.1). Also, 
the substrate emittance is constant (h { =h u =\), no reflectance at inter- 
faces ( Pv0 = p vs = 0), and &E/kT E =0.\. (a) No scattering 

(n y =n w =n g =0); (b) equal scattering and absorption coefficients 

(n / =n tf =n g =o.5). 

5(b)] are more than a factor of 2 less than those for 
f u = f 0.01 [Fig. 5(a)]. As stated earlier, the rare-earth 
oxides are expected to have f u and // values between 0.01 
and 0.1. 

The dependence of efficiency on optical depth for the 
case of equal scattering and absorption coefficients 
(n,= n„=n^=0.5) is shown in Fig. 6. Comparing Figs. 5 
and 6 shows that scattering has negligible effect on the 
magnitude of y E . Scattering merely results in the maxi- 
mum y E occurring at larger K g than for the nonscattering 
case. These results were also deduced from Figs. 3 and 4. 


Now consider the importance of the interface reflectiv- 
ities p x<} and p vs . Equation (48) shows that the major con- 
tribution to e v of the vacuum-film interface reflectivity is 
the ( 1 — p v0 ) term. Therefore, if p x<) is a constant this term 
will cancel in the ratios e/e g and € u /e g that appear in y E 
[Eq. (55a)]. As a result, the p vs term in Eq. (48) will have 
the major effect on r\ E rather than p x<} . 

To determine just how important p ^ and p vs are in 
determining rf E we approximated the reflectivities in the 
following manner. For the vacuum-film interface the nor- 
mal reflectivity for a specular reflective dielectric interface 
was used: 17 



Also, since n f > 1 the maximum incidence angle d M that 
allows radiation to escape the film is no longer u/2 but is 
given by Eq. (13). There are limited data available on rij- as 
a function of wavelength for rare-earth oxides. Marcinow 
and Truszkowska 20 found that n f for Gd 2 0 3 , Eu 2 0 3 , 
Er 2 0 3 , and Yb 2 0 3 is nearly constant for 0.5<A<2.5 pm. 
Also, they found that 1.85</iy<1.95 for these four rare- 
earth oxides in the 0.5<A<2.5 pm wavelength region. The 
index of refraction data of Refs. 21 and 22 are in agreement 
with these results. Therefore, for the calculations a con- 
stant value of rif = 1.9 was chosen. As a result, p v0 = 0.096 
and Pm— 0.85. As already discussed, the substrate must be 
a low-emittance (e^<0.05) material. In a later discussion of 
substrate materials rhodium and copper are recommended 
as substrate materials. Since both copper and rhodium are 
opaque, p vs +e vs = 1; therefore, p vs — 1 — e vs was used in the 
calculations. 

In Fig. 7 the effect of p ^ and p vs on the optimum 
efficiency is presented for the same conditions as Fig. 3. 
Figure 7(a) is for no scattering and Fig. 7(b) is for equal 
scattering and absorption coefficients. Comparing Figs. 3 
and 7 shows that there is about a 10% increase in y E for 

f opt 

e^=0.01, but at £^=0.17, tj e has increased by nearly a 
factor of 2 for n f= 1.9. Also note that the optimum optical 
depth is approximately a factor of 3 smaller. As discussed 
earlier this change results from p vs rather than p ^ for a 
constant rif.Xirif were strongly dependent on wavelength 
then the (1— p v ) term in the emittance [Eq. (48)] would 
also depend on wavelength. Thus, the emittance ratios 
e/€ g and e u /e g that appear in y E [Eq. (55a)] would also 
depend on wavelength. Therefore, in that case y E could be 
significantly effected by p v0 . 

IV. EMISSIVE POWER OF THIN-FILM SELECTIVE 
EMITTER 

Other than efficiency, the power emitted in the emis- 
sion band is the most important performance parameter for 
a thin-film selective emitter. This emitted power/area P E is 
given by the numerator of Eq. (55). Now define the di- 
mensionless power density as follows: 

P E 6 8 / E - t±EJ2 e Eb dE 

Pe=—^= dF ’ (62) 

°SB *E J0 e Eb“E 
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TABLE I. Emission-band data for rare-earth oxides. 


Emitter material 

Emission-band 

transition 

Photon energy 
at center of 
emission band 
E g (eV) 

Photon wavelength 
at center of 
emission band 

(jum) 

Dimensionless bandwidth 

Emitter temperature 
for maximum efficiency, 

T e/ opt =4 eEg/fc 

<K) 

Ytterbia 

Yb 2 0 3 

2 F 5/2 to 2 F m 

1.29 

0.955 

0.18 

3740 

Erbia 

Er 2 0 3 

4 1 11/2 t0 4 7 15/2 

0.827 

1.5 

0.05 

2400 

Holmia 

Ho 2 0 3 

5 / 7 to 5 / 8 

0.62 

2.0 

0.10 

1800 

Neodymia 

Nd 2 0 3 

4 A 1/2 4 ^9/2 

0.52 

2.4 

0.15 

1500 


where cr SB is the Stefan-Boltzmann constant 
(or SB = 5.67X 10 -12 W/cm 2 K 4 ) and 7^ is the emitter tem- 
perature. Using Eq. (6) for the blackbody emissive power 


e ej, gives the following: 


15 p(l + f/2) x 3 

Pe—~J i dx > 

it Js(\-t/i) er—i 

(63) 

where 


E g A Eg 

c — t~ * 

kT E ' E g 

(64) 


The power density p E behaves differently, as a function of 
optical depth, than the efficiency. As Figs. 5 and 6 show, 
there is an optimum optical depth for maximum efficiency. 
However, as Eq. (63) shows, p E depends on optical depth 
only through e g . And as Fig. 2 shows e g is a monotonically 
increasing function of optical depth. As a result, p E will 
also be a monotonically increasing function of optical 
depth. Also, the substrate emittance e sg will affect p E only 
for small optical depths. Most important, however, is the 
effect of scattering on p E . As discussed earlier, scattering 
has a negligible effect on the optimum efficiency; however, 
scattering reduces p E since e g decreases with increasing Ct g 
(Fig. 2). These points are illustrated in Fig. 8 where p E is 
shown as a function of optical depth. Figure 8(a) is for the 
nonscattering case (fl g =0) and Fig. 8(b) is for equal scat- 
tering and absorption coefficients (H g =0.5). Notice that 
p E increases more rapidly for the nonscattering case. Near 
maximum p E is attained at K g > 1; however, with scattering, 
AT g >2 for maximum p E . 

Figure 9 shows p E as a function of Eg/kT E for the same 
conditions as Fig. 8 with K g = 1 for the nonscattering case 
[Fig. 9(a)] and K g = 2 for ft g =0.5 [Fig. 9(b)]. With these 
values of K g the power densities will be close to their max- 
imum values. Similar to rj E , p E has its maximum value at 
Eg/kT E ~ 4; however, the location of maximum r) E and p E 
will shift if /, g, and h are different below and above the 
emission band. This is illustrated in Fig. 10 where v E and 
p E are shown as functions of Eg/kT E for f u = 0.01, f { —0A, 
g B =g/=A B =Ai = l, &Eg/E g =0A with no scattering 
(fl g =0) or reflectance at the interface and for two values 
of substance emittance {e sg — 0.01, 0.05). As stated earlier, 
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when f f/^f u the optimum optical depth is a function of 
Eg/kT E . Therefore, each point on the curves in Fig. 10 
corresponds to a different optical depth. The power density 
p E was calculated using the optimum optical depth. As can 
be seen, if / )> f u the maximum efficiency and power den- 
sity are shifted to a lower value of Eg/kT E than is the case 
where f u = // (Figs. 3 and 4); however, the dependence of 
maximum p E on Eg/kT E is much less than the maximum 
efficiency dependence. Also note that p E is significantly 
effected by the substrate emittance compared to Fig. 9 
where K g is a constant. 

V. POSSIBLE SUBSTRATES FOR RARE-EARTH 
OXIDE THIN-FILM SELECTIVE EMITTERS 

From the analytical results just discussed we know 
that a thin-film selective emitter must have a low-emittance 
substrate in order to have high efficiency; also, the sub- 
strate must have a high melting temperature. These two 
requirements greatly limit the choice for a substrate mate- 
rial. Copper, silver, and gold all have low emittance 23 
(e s <0A) y however, their melting points are also low. Cop- 
per with a melting point of 1356 K may be suitable for use 
with Nd 2 0 3 . For the high temperatures required by Er 2 0 3 , 
Ho 2 0 3 , and Yb 2 0 3 two possible substrates are sapphire 
with a melting point of 2326 K and rhodium (Rh) with a 
melting point of 2233 K. Sapphire has low emittance 
(e vs <0.1) for wavelengths A <4.0 pm, however, for 
5<A<15 pm the emittance is large. 23,24 Therefore, sapphire 
is most suitable for the highest temperatures (7>2000 K) 
where most of the radiation is at A < 4 pm. Rhodium has 
total emittance £ 5 <0.1 up to the highest temperature 
(1500 K) for which data are available. 23 As a result rhod- 
ium should be a suitable substrate for all the rare-earth 
oxides. 

VI. CONCLUSION 

Using a radiation transfer analysis that includes isotro- 
pic scattering an expression for the spectral emittance of a 
thin film was developed. Interference effects were ne- 
glected. Splitting the spectral optical properties of the rare- 
earth oxides into three regions (below emission band, emis- 
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FIG. 5. Dependence of emitter efficiency rj E on emission-band optical 
depth A^for no scattering (ft g =n u = fl/=0) and constant substrate emit- 
tance ( h { —h u = 1 ). Also, E/kT £ =4, &Eg/E g =OA, and no reflectance at 
interfaces (p^ = p vs —0). (a) Both below- and above-emission-band op- 
tical depths equal to 0.01 (emission-band optical depth) (f u — //= 0.01 ); 
(b) both below- and above-emission-band optical depths equal to 0.1 
(emission-band optical depth) (/,=/„ =0.1). 


sion band, and above emission band) and using the 
spectral emittance expression, the emitter efficiency rj E was 
calculated. Calculated optimum efficiencies range from 
t] e ~0.35 to 7] e ~0.1 . 

Two important results of the analysis are the following. 
First, in order to achieve high efficiency the film substrate 
must have very low emittance (e 5 <0.05). Second, there is 
an optimum optical depth to obtain maximum rj E . For the 
scattering coefficient cr v less than the or equal to the ab- 
sorption coefficient a vf scattering has a negligible effect on 
the optimum efficiency. However, the optical depth re- 
quired to obtain maximum efficiency with scattering is 


FIG. 6. Dependence of emitter efficiency t] e on emission-band optical 
depth K g for equal scattering and absorption coefficients 
(n^=n u = n/=0.5) and constant substrate emittance {h u ~h ( = 1). Also, 
&Eg/E g =0A and no reflectance at interfaces (p v0 = p vs =0). (a) Both 
below- and above-emission-band optical depths equal to 0.01 (emission- 
band optical depth) (/ tt =/,=0.01 ); (b) both below- and above- 
emission-band optical depths equal to 0.1 (emission-band optical depth) 
(//=A=o.i). 


greater than the optical depth required to obtain maximum 
efficiency without scattering. Also, scattering reduces the 
emissive power density. Finally, the maximum efficiency 
occurs when Eg/kT where E g is the photon energy at 
the center of the emission band and T E is the emitter tem- 
perature. 

The low-emittance and high-temperature requirements 
on the substrate severely limit the choice of substrate ma- 
terial. Copper may be suitable for use with Nd 2 0 3 . For the 
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Eg/XTE 



V kT E 



Kg 



FIG. 8. Dependence of dimensionless emitter power density , 
pE=^£/ a sh^Et on emission-band optical depth K g for substrate emittance 
e 5Jf =0.01 and 0.05. Also, E^kT E — 4, LEg/E g =0.\ t and no reflectance at 
interfaces (p l o=p«=0)- ( a ) No scattering case (£l g ~0) and K g =\.0; 
(b) equal scattering and absorption coefficients (ft^0.5) and K g =2.0. 


, e vb n v ( q + ( 0) 

S v (K v )=nl(l-n v ) —+—[E 2 (K v ) 

J 7 T Z \ 7 T 

q-( 0) rv 

+£ 2 (tf w -tf v ) Hh 5 V (^J) 

^ Jo 

X£,(|A*-A V | )<«*). (Al) 


FIG. 7. Effect of finite reflectance at vacuum-film interface and film- 
substrate interface on optimum efficiency VE op{ - Pvo given by the re- 
flectance for normal incidence [Eq. (61)], n f = 1.9, p^—0.096, and 
p vs — 1 — € sg . Also, f i= f u — 0.01, gi=g u = 1, hi~h u — l, and 1. 

(a) No scattering (ft^ft^ft/^0); (b) equal scattering and absorption 
coefficients (D u — n # =n/=0.5). 


other rare-earth oxides, sapphire and rhodium are possible 
substrate choices. 


APPENDIX: SOURCE FUNCTION SOLUTION 

Equation (3) must be solved for the source function 
S V (K V ) that appears in the radiative transfer Eqs. ( 1 ) and 

(2) . Since the boundaries are assumed to behave diffusely, 
/+(0,/z) and i~(K ve ,—fi) are independent of fi and Eq. 

(3) becomes the following: 


In order to reduce this integral equation for S v to a differ- 
ential equation, the exponential integrals are approximated 
as follows: 

E n (x) =a n exp( — b„x). ( A2) 

For E x {x) the values a I =Z? 1 = 2 were used in calculations, 
which is the same approximation used in Ref. 18. It should 
be noted that E x (0)^oo; however, the approximation 
given by Eq. (A2) for/s^x) yields ^(O) =a ] . For E 2 (x) 
the approximations a 2 = 3/4, b 2 = 3/2 were used. 17 

To convert Eq. ( Al ) to a differential equation, use the 
approximations given by Eq. (A2) for E x and E 2 and then 
take the second derivative with respect to K v . This yields 

d 2 S v 

— 2— a,S v = -a 2 +a 3 ex P( ~b 2 K v ) +a 4 exp(M^v). 
aK v 

(A3) 

where 
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FIG. 9. Dependence of dimensionless emitter power density, 
p E = on dimensionless emission-band energy Eg/kT E for sub- 

strate emittance e sg =0.01 and 0.05. Also, AEg/E g — 0. 1 and no reflectance 
at interfaces (p^^p^—0). (a) No scattering (tl g ~ 0) and K g = 1.0; (b) 
equal scattering and absorption coefficients (0^=0. 5) and K g =2.0. 



Eg/KTE 



E^T e 


FIG. 10. Effect of unequal below-emission-band optical depth f t and 
above-emission-band optical depth /„ (//=0.1,/ u = 0.01 ) on optimum 



(A4) 

efficiency r] E and dimensionless power density p E . Also, no scattering 
(n u = n/=ft^=0), constant substrate emittance {h u =h/= 1), 
&Eg/E g =Q.\ t and no reflectance at interfaces (p^ 0 = p Vf — 0). 

a 2 = «}(l-O v )^\ 

J 77 

( A5) 


«3-- 2 (*,-* 2 ) v > 

(A6) 

c, = ( 0 ) -^ 29v - (*«,) +^ 3 «^ vi ]e- , 1 

(A9) 

a 2^v / l 2 El\ $ (**> / l is \ 

« 4 -- ^ (bi—b 2 ) exp {-b 2 K vd ). 

L 77 

(A7) 

C 2 =^ [AjffUO )-A x q;{K vd )+A i n 1 J e vb \, (A10) 


The solution to this linear second-order differential equa- 
tion is 

S V (K V ) = C{ exp( yjcciK v ) + C 2 exp( — 

+ C 3 exp ( - b 2 K v ) 4- C 4 exp (b 2 K v ) + C 5 , 

(A8) 

where the constants (C { , C 2 , C 3 , C 4 , and C 5 ) must satisfy 
Eqs. (Al) and (A3). If Eq. (A8) is substituted in Eq. 
(A3) then results for C 3 , C 4 , and C 5 are obtained. Using 
those results for C 3 , C 4 , and C 5 and then substituting Eq. 
(A8) in Eq. (Al) will yield (after a great deal of algebra) 
results for C x and C 2 \ 


«3 

^ v +(0) , 

(All) 

b\-a x 

O- 

II 

X e -y q -{K vd ), 

77 

(A12) 

1 i-n v 

n 2 je vb . 

(A13) 

7T 1 — Z 


The quantities A } , A 2f A 3 , A 4 , r, y , z, and y are given by 
Eqs. (32)-(41). For no scattering, Cl v —z=y— A a — 0 and 
all the Cs vanish except C 5 . Thus, as expected, only the 
blackbody emissive power e vb contributes to the source 
term. 
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